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I. INTRODUCTION 
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We develop the real stabilization method within the framework of the reiativistic mean field 
(RMF) model. With the self-consistent nuclear potentials from the RMF model, the real stabilization 
method is used to study single-particle resonant states in spherical nuclei. As examples, the energies, 
widths and wave functions of low- lying neutron resonant states in 120 Sn are obtained. These results 
are compared with those from the scattering phase shift method and the analytic continuation in 
the coupling constant approach and satisfactory agreements are found. 
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0^ ' The investigation of continuum and resonant states is an important subject in quantum physics. In recent years, 
there has been an increasing interest in the exploration of nuclear single particle states in the continuum. The 
construction of the radioactive ion beam facilities makes it possible to study exotic nuclei with unusual N/Z ratios. 
In these nuclei, the Fermi surface is usually close to the particle continuum, thus the contribution of the continuum 
and/or resonances being essential for exotic nuclear phenomena 0, 0, H, EL B 0- It has been also revealed that the 
contribution of the continuum to the giant resonances mainly comes from single-particle resonant states 0, H| ■ 

For the theoretical determination of resonant parameters (the energy and the width), several bound-state- like 
methods have been developed. The complex scaling method (CSM) describes the discrete bound and resonant states 
on the same footing Q. In this method, a complex coordinate scaling is introduced to rotate the continuum into 
the complex energy plane and the wave functions of resonant states, but not scattering states, are transformed into 
square-integrable functions [lfj| | . Although it involves the solution of a complex eigenvalue problem which causes some 
difficulties in practice, the CSM has been widely and successfully used to study resonances in atomic and molecular 
systems [ll|, O HI] and atomic nuclei 0, [13, EH EH E3| ■ The analytical continuation in the coupling constant (ACCC) 
approach is based on an intuitive idea that a resonant state can be lowered to be bound when the potential becomes 
more attractive or cquivalcntly the coupling constant stronger, thus a resonant state being related to a series of 
. bound states via an analytical continuation in the coupling constant [l?], E3, E3] ■ Combined with the cluster model, 
the ACCC approach has been used to calculate the resonant energies and widths in some light nuclei 0, [2l[ . An 
attempt to ex plo re the unbound states by the ACCC approach within reiativistic mean field (RMF) model was first 
■ made in Ref. [22J where resonant parameters of some low lying resonant states in 16 O and 48 Ca obtained from the 
ACCC calculations are comparable with available data. The wave functions of nuclear resonant states were also 
determined by the ACCC method where the bound states are obtained by solving either the Schrodinger equation 



with a Woods-Saxon potential [23( or the Dirac equation with self-consistent RMF potentials [24 1 



The real stabilization method (RSM) is another bound-state- like method [25| . The equation of motion of the system 
in question is solved in a basis (2f| or a box (26| of finite sizes, thus a bound state problem being always imposed. The 
RSM uses the fact that the energy of a "resonant" state is stable against changes of the sizes of the basis or the box. 
It has been used to calculate the resonance parameters in elastic and inelastic scattering processes 127 1 12811 . Some 
efforts have also been made in order to calculate more efficiently resonance parameters with the RSM [23, l3C)L l3ll |32| . 
In this work, we investigate single particle resonances in atomic nuclei by combining the RSM and the reiativistic 
mean field (RMF) model. 

The paper is organized as follows. In Sec. [Til we give briefly the formalism for the RSM and the RMF model. The 
numerical details, the results for 120 Sn and discussions arc given in Sec. IIIII Finally we give a brief summary. 
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II. FORMALISM OF THE RMF MODEL AND THE RSM 



A. The relativistic mean field model 



The basic ansatz of the relativistic mean field (RMF) model is a Lagrangian density where nucleons are described 
as Dirac spinors which interact via the exchange of several mesons (a, u>, and p) and the photon [33l. |34|. l35l. l36l. l37j| . 



ipi (iP - M)ipi + -d^ad^a -U(a)- g a ipiVipi 



-R flv R' J,v + ^m 2 p p p pf - gptpipri/ji 



(1) 



Tx^ = j^x*, M the 



where the summation convention is used and the summation over i runs over all nucleons, = = 
nuclcon mass, and m a , g a , m u , g u , m p , g p masses and coupling constants of the respective mesons. The nonlinear 
self-coupling for the scalar mesons is given by [38[ 
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and field tensors for the vector mesons and the photon fields are defined as 

On. 

dp,p v - dvp^ - g p (p^ x p v ), 



(2) 



(3) 



F„ 



d^Av — d v A p . 



The classical variation principle gives equations of motion for the nucleon, mesons and the photon. As in many 
applications, we study the ground state properties of nuclei with time reversal symmetry, thus the nucleon spinors 
are the eigenvectors of the stationary Dirac equation 



[a-p + V(r) + P(M + S(r))} ^(r) = e^(r), 

and equations of motion for mesons and the photon are 

" (-A + d a U(a))a(r) = -g aPs (r), 

(-A + ml)ui°{r) = g u p v (r), 

(-A + m 2 p ) p°(r) = g p p 3 {r) 7 

—AA°(r) = ep p (r), 



(4) 



(5) 



where ui° and A are time-like components of the vector oj and the photon fields and p° the 3-component of the 
time-like component of the iso-vector vector p meson. Equations dU) and ([5]) are coupled to each other by the vector 
and scalar potentials 



V(r) = g^°(r)+g p T 3 p (r) + e 
S(r) = g a cr{r), 



1-T 3 



A°(r), 



(6) 



and various densities 



pv{r) = EiLiV'iM^iO"), 



(7) 



Pc{r) = Ei=iV>lM 
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For spherical nuclei, meson fields and densities depend only on the radial coordinate r, the Dirac spinor reads 



jm\ 



F K (r) 



(8) 
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with Y 



the spin spherical harmonics. The radial equation of the Dirac spinor, Eq. ((4]), is reduced as 



e aGa — 
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The meson field equations become simply radial Laplace equations of the form 

-)2 



are the meson masses for 



= (j,lo, p and zero for the photon. The source terms are 
s<p(r) = 



-g<rps(r) - g2<J 2 {r) - g 3 cr 3 {r), for a, 

g u p v {r), for u, 

9 P P3(r), for p, 

> ep c (r), for A, 



with 



47 rr 2 p s (r) = £f =1 (|Gi(r)| 2 - |W| 2 ), 
4^(0 = Eti(I^W| 2 + l^(r)| 2 ), 



47rr 2 /r?3(r) 
47rr 2 p c (r) 



Eti2^(|G,(r)| 2 + |^(r)| 2 ), 
Ef=i(|-*i)(l^(r)| 2 + |i^(r)| 



(9) 



(10) 



(11) 



The above coupled equations can be solved iteratively in r space [39( or in the harmonic oscillator basis [4C 
no sea and the mean field approximations. 



(12) 



using the 



B. The real stabilization method in coordinate space 

With the self consistent vector and scalar potentials V(r) and S(r), the Dirac equation ([9]) is solved in a spherical 
box of the size i? m ax under the box boundary condition, and thus the continuum is discrctizcd. When i? m ax is large 
enough, the energy of a bound state does not change with i? max . In the continuum region, there are some states 
stable against the size of the box, i.e., the energy of each of such states is almost constant with changing i? max ; such 
stable states correspond to resonances. 

The resonant parameters, E 1 and T, may be obtained by fitting the energy E and the phase shift T)(E) in an energy 
range around a resonance to the following formula [25| . 

Vi(E) = m, P ot(E) + tan- 1 ( -^-) ■ (13) 



E — E-y 



The phase shift rji (E) can be calculated as |41| 



r Rmax Xi{r) [E - H(r)] f(r) sinkrdr 

tan [rji —J = ^ , (14) 

Jo ™ x Xi(r) [E - H(r)] f(r) cos krdr 
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with f(r) satisfying /(r) — > 1 when r — > 00 and /(0) = /'(0) = 0. However, Eq. ([M]) converges very slowly with the 
box size due to the influence of the non zero centrifugal potential at large fi m M [4l|. 

In the present work, we use a simpler method proposed by Maier et al. [26( in which it's not necessary to calculate 
the phase shift. The resonance energy is determined by the condition d 2 E/dB% aax = and the corresponding box 
size is labeled as R max , i.e., E y = E(R max ). The width is evaluated from the stability behavior of the positive energy 
state against the box size around i? max - 

When r is large enough, the nuclear potentials S(r) and V(r) vanish, G(r)/r satisfy 

d 2 G , ( 2 k(k+1) 



dr 2 



a 2 - v „ ; G = 0, (15) 



with a 2 = E 2 — M 2 . The general solution reads 

G{r) oc ar [cos r/i ji(ar) — sin 77; ni(ar)] , (16) 
When r — > oo, G(r) oc sin(or — ^ + 77/). Therefore when i? max is large enough, 

In 

a-Rmax - y + Vl = n^. (17) 

Under the assumption that the phase shift from the potential scattering r]^ po t(E) varies slowly with respect to the 
box size, i.e, dr]i_ po t/dR max ~ 0, one derives from Eqs. (|T3|) and (fT7|) the formula, 



2JE 2 + 2E^M 



T= i -r. (18) 

-(B 7 + M)i? roax - (£ 2 + 2 £ 7 M) [a^/Sii^l^J 



In the non-relativistic limit, E y -C M, Eq. (fT8|) is reduced to 



2 V / ^ 7 (19) 



— -Rmax ^ 2i? 7 [d^/di?max|fl max ] 

which is essentially the same as Eq. (8) in Ref. [26[ except for that here natural units with h = c = 1 is used. 

III. RESULTS AND DISSCUSION 

In this section we present the results of the RSM in the framework of the RMF model. In our calculation, we use 
for the Lagrangian density the effective interactions PK1 [12] and NL3 (43J. We change the size of the box in a large 
range (7 fm < i? m ax < 60 fm if not specified) in order to find not only narrow resonances but also wide ones. We 
take 120 Sn as an example and compare the results for neutron resonances from the RSM with those from the ACCC 
approach [24| and the scattering phase shift method [44| . 

By examining the stability of low lying positive energy states against the size of the box i? ma x, we find that except 
for s state, there are neutron resonances in 120 Sn with the orbital angular momentum I up to 6. 

The narrowest resonance is an i 13 / 2 state which lies at about 3.45 MeV above the threshold. The positive energy 
VI13/2 states in a box of different sizes are shown in Fig. [T] With the box size i? max increasing, the lowest ^ii 3 / 2 state 
first falls down quickly then its energy becomes constant in a large region of i? m ax- After it crosses with the second 
lowest z/ii3/2 state at around i? m ax = 25 fm, the lowest one falls down again. Similar level-crossings occur regularly 
at larger i? ma x- This stability behavior implies that there is a narrow 1^13/2 resonant state with approximate energy 
3.45 MeV. The resonant energy 3.469 MeV and -R max = 20.1 fm are obtained under the condition d 2 E/dR 2 aax = 




FIG. 1: Positive energy ^13/2 states in Sn under different box boundary conditions. 
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FIG. 2: The radial wave function for the neutron id 13 / 2 resonant state in 120 Sn. 



TABLE I: Energies and widths of the single neutron resonant state 113/2 in 120 Sn from different E ~ i? ma x curves (cf. Fig. [TJ. 
fl max is in fm and E 1 and T are in MeV. 



Rmax 




100 xT 


Rmax 




100 xT 


20.1 


3.4688 


0.327 


47.2 


3.4686 


0.470 


30.3 


3.4686 


0.426 


55.3 


3.4685 


0.478 


38.9 


3.4686 


0.456 


63.2 


3.4686 


0.483 



from the first E ~ i? max curve (labeled as "lii3/2") in Fig.Q] The width 0.003 MeV are obtained from Eq. (fT5|) . One 
of the approximations made in deriving Eq. (jT5J) is i? max should be large. We next examine the dependence of the 
resonant parameters on the box size by calculating E 1 and T from other E ~ _R max curves with larger i? max . For 
this purpose the calculations with -R max up to 65 fm are carried out. The results are given in Table [TJ The energy is 
almost a constant with increasing R max . The variation between the widths obtained from adjacent E ~ i? max curves 
decreases with _R max and is about 1% at -R max ~ 63 fm. For other resonant states presented in this work, we also 
make similar investigations. Once it converges to within 1%, the value of the width is assigned to a resonant state. 
The wave function of the resonant state v'113/2 is given in Fig. [2] From this figure one can find that this state is almost 
localized inside the nucleus which is consistent with the small width T. 

Although it lies below the state ^13/2, the resonant state ^5/2 is about an order of magnitude wider than zdi3/2- 




FIG. 3: Positive energy i^f 5 / 2 states in 120 Sn under different box boundary conditions. 
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FIG. 4: The radial wave function for the neutron vi*,/2 resonant state in Sn 
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FIG. 5: Positive energy fin/2 states in Sn under different box boundary conditions. 



The reason is that the centrifugal barrier for ffs/2 (I = 3) is much lower than that for 1^13/2 (2 = 6). The stability 
behavior of the state ^5/2 is presented in Fig. [31 The resonant energy and width are 0.870 MeV and 0.064 McV, 
respectively. The wave function of this state is also well localized as shown in Fig. 2J 

Except the state fii.3/2, there is another neutron resonant state with I = 6, i.e., the fin/2 state. As shown in Fig. [5l 
there is a less stable state lying at about 10 McV. Although it shares the same centrifugal barrier with ^13/2, fin/2 
is much wider because its energy is much larger than that of fii3/2- In Fig. [6j one finds also that the wave function 
of M13/2 oscillates very much even at r = 50 fm. The resonance parameters for this state are E 1 = 9.811 MeV and 
r = 1.275 MeV respectively. 

The energies and widths of single particle neutron resonant states obtained from the RSM calculations are sum- 
marized in Table [ill We also calculate these resonances using the parameter set NL3 for the Lagrangian density in 
the RMF model and compare the present results with those from the ACCC approach and the scattering phase shift 
method [45| in Table HT1 In Fig. [3 the comparison is also made in a planar E 1 -T plot. For the two low lying resonant 
states, vt§/2 and ^13/2, the three methods give consistent results both for the energy and the width. Although similar 
energies are obtained from these three models for higher resonances, fin/2 an d fjis/2; clear differences occur among 
the widths from the ACCC approach and the scattering method and the results from the RSM lies in between. 
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FIG. 6: The radial wave function for the neutron vin/2 resonant state in 120 Sn. 
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FIG. 7: Energies and widths of single neutron resonant states in Sn from different methods. RMF-ACCC, RMF-RSM, and 
RMF-S represents results from the analytical continuation in the coupling constant approach, the real stabilization method, 
and the scattering phase shift method in the framework of the RMF model with NL3 parameter set. 



TABLE II: Energies and widths of single neutron resonant states in °Sn from different methods. RMF-ACCC, RMF-RSM, 
and RMF-S represents results from the analytical continuation in the coupling constant approach, the real stabilization method, 
and the scattering phase shift method in the framework of the RMF model. In the RMF-RSM calculations, both PK1 and NL3 
parameter sets are used. All quantities are in MeV. 





RMF-RSM (PK1) 

e 1 r 


RMF-RSM (NL3) 

E-y r 


RMF-ACCC (NL3) 

e 1 r 


RMF-S (NL3) 

e 1 r 


^13/2 
fill/2 
fjl5/2 


0.870 0.064 
3.469 0.005 
9.811 1.275 
12.865 1.027 


0.674 0.030 
3.266 0.004 
9.559 1.205 
12.564 0.973 


0.685 0.023 
3.262 0.004 
9.60 1.11 
12.60 0.90 


0.688 0.032 
3.416 0.005 
10.01 1.42 
12.97 1.10 
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IV. SUMMARY 



In summary, the real stabilization method (RSM) has been developed within the framework of the relativistic mean 
field (RMF) model. With the self-consistent nuclear potentials provided by the RMF calculations with the parameter 
sets PK1 and NL3 for the Lagrangian density, the Dirac equation for the neutron is solved in the coordinate space 
under the box boundary condition. By investigating the stable behavior of the positive energy states against changes 
of the box size, the resonant states are singled out. The RMF-RSM is used to study single neutron resonant states in 
spherical nuclei. As examples, the energies, widths and wave functions of low-lying neutron resonant states in 120 Sn 
are obtained. Since a very large box size is used, even wider resonances can also be found. These results are compared 
with those from the scattering phase shift method and the analytic continuation in the coupling constant method and 
satisfactory agreements are found. 
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